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Abstract 

We theoretically study the trapping time distribution and the efficiency of the excitation energy 
transport in dendritic systems. Trapping of excitations, created at the periphery of the dendrimer, 
on a trap located at its core, is used as a probe of the efficiency of the energy transport across 
the dendrimer. The transport process is treated as incoherent hopping of excitations between 
nearest-neighbor dendrimer units and is described using a rate equation. We account for radiative 
and non-radiative decay of the excitations while diffusing across the dendrimer. We derive exact 
expressions for the Laplace transform of the trapping time distribution and the efficiency of trapping 
and analyze those for various realizations of the energy bias, number of dendrimer generations, 
and relative rates for decay and hopping. We show that the essential parameter that governs the 
trapping efficiency, is the product of the on-site excitation decay rate and the trapping time (mean 
first passage time) in the absence of decay. 
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I. INTRODUCTION 



During the past decade dendritic molecular systems or dendrimers have received consid- 



erable attention 



Dendrimers are synthetic highly branched tree-like macro- 
molecules consisting of a core and several branches which are not connected geometrically 
to each other and are built self-similarly. Theoretically the process of building a dendrimer 
can be repeated ad infinitum to obtain a dendrimer with any number of generations, but 
in practice this number is currently limited to fifteen [6j. The first three dendrimers with 
coordination number z = 3 are schematically depicted in Fig. ^ 

Dendrimers hold great promise for creating artificial light harvesting systems. Indeed, 
because of their branched nature, the number of units at their periphery grows exponentially 
with the number of generations . Therefore, if light absorbing states are located at the pe- 
riphery, the cross-section of absorption grows exponentially with the number of generations. 
Combining this with a possiblvefficient energy transport to the core may result in efficient 
light harvesting systems 

A growing flux of publications exists, bot 



16, h/Zll and theoretical 



21 




reporting on op- 



tical and transport properties of dendritic systems. For instance, in Ref. |8| compact and 
extended polyphenylacetylene dendrimers in solution were studied experimentally. Inter- 
preting their spectroscopic measurements, the authors argued that the optical excitations 
in these systems are localized on dendrimer subunits. They also made the important point 
that extended dendrimers were characterized by a so-called energy funnel: the excitation 
energies of dendrimer units decreases from the periphery towards the core, thus providing an 
energetic bias for transport of the excitation towards the core. Such a bias is favorable for 
efficient transport of the absorbed energy across the dendrimer. Using quantum chemical 
calculations, the existence of an energy funnel towards the core was further confirmed in 
Ref . I2I . Experimental investigations of the energy transport in these systems Q revealed 
the fact that it occurred through a multi-step (incoherent) hopping process with an effi- 
ciency of 96%. On the other hand, low-temperature measurements of the energy transport 
in distyrylbenzene-stilbene dendrimers with a nitrogen core gave clear indications of coher- 
ent interactions between the dendrimer subunits ilTil. In this class of organic dendrimers, 

n 

ultra-fast higher-order nonlinearities were also reported which makes them potentially 
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promising systems for applications in nonlinear optical switching elements. 

The first theoretical efforts on dendrimers were mostly focused on analyzing the mean 
first passage time, i.e., the average time it takes for an excitation, created somewhere at 
the periphery of the dendrimer, to reach its center, where a trap is located. This problem 
was addressed in much detail in Refs. 3, IS, where the mean first passage time was 
calculated both in the presence and the absence of a fixed energy bias. The presence of a 
random bias was considered in Ref. js^, where it was found that the randomness tends to 
reduce the transport efficiency. Large-scale numerical calculations of the statistics of the 
number of sites visited by the excitation and its mean square displacement in very large 
dendritic structures were performed in Ref. 23 1, while the authors of Ref. discussed the 
kinetics of symmetric random walks in compact and extended dendrimers of a small number 
of generations (< 4). 

In all papers cited above, the basis of the analysis was the assumption of incoherent hop- 
ping transport between dendrimer units. However, the actual nature of the excitation trans- 
port and optical dynamics is still under debate. Thus, the alternative approach of coherent 
excitons extending over several dendrimer units has also been considered, in particular to de- 
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and the enhancement 



scribe energy transport from the periphery to the core 
of the third-order optical susceptibility 29]. An other interesting branch of experimental and 
theoretical studies of dendritic systems concerns the dynamics of dendrimer-based networks, 
i.e., extended systems using dendrimers as building blocks |3flL l3lL Recently, such 

networks have attracted much attention due to their two levels of structural organization. 

Among the factors that govern the efficiency of energy transport in dendritic systems, 
such as the energy bias and the presence of disorder, the radiative and non-radiative decay 
of optical excitations during their random walk to the core represents an additional channel, 
especially for dendrimers with a large number of generations. To the best of our knowledge, 
this factor has thus far not been addressed in the literature. Moreover, the statistics of the 
trapping time has not been studied in detail either. In this paper, we intend to fill these 
gaps and show that intimate relations exist between the analysis of both issues. We will 
show that in the presence of decay, the key parameter that governs the trapping efficiency 
is the product of the on-site excitation decay rate and the mean first passage time in the 
absence of decay. 

The outline of this paper is as follows. In Section |n] we present our model, which is based 
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on incoherent motion of excitations across the dendrimer, described by a rate equation that 
accounts for both trapping at the core and excitation decay during the random walk towards 
the core. We also discuss several quantities relevant to the problem under study, such as 
the distribution of survival times, the mean survival time, and the efficiency of trapping. 
Section ITTTl deals with deriving exact expressions for these quantities in the Laplace domain. 
In Sec. IIV| we present a detailed analysis of the trapping time distribution in the limit of 
vanishing excitation decay. Effects of the excitation decay on the trapping efficiency under 
various conditions (sign of the energy bias, number of generations, ratio of decay and hopping 
rates) are discussed in Sec. El Finally, we conclude in Sec. IVII 

II. MODEL AND PERTINENT QUANTITIES 

As was already mentioned in the Introduction, we treat both the motion of the excitation 
over the dendrimer units and the trapping at the core as incoherent nearest-neighbor hopping 
processes. We label the dendrimer units (sites) by the index i {1 < i < A/"), where A/" is the 
number of units, while i = denotes the trap located at the core (cf. Fig. 1). The trap is 
considered irreversible, i.e., once the excitation hops onto it, it never returns to the body of 
dendrimer. Then, the system of rate equations for the excitation probabilities (populations) 
Pi of the dendrimer units reads 



Here, the dot denotes the time derivative, the summation performed over sites j 

that are nearest-neighbors of the site i, the prime in the second summation of Eq. ()lb|) 
indicates that j 7^ 0, 7 is the exciton decay rate (assumed independent of i), and kij is the 
rate of hopping from site j to site i, including i = 0. The hopping rates meet the principle 
of detailed balance: kij = kjiexp[{Ej — Ei)/kBT], where Ei is the excitation energy of the 
site i, ks is the Boltzmann constant, and T is the temperature. In this sense the energy of 
the trap, Eq, is considered infinitely low. Initially, the excitation is outside the trap, i.e., 
]9o(0) = 0, while one of dendrimer units is excited, Pi{0) = da^. 

The quantity r{t) = po represents the instantaneous trapping rate and will be used to 




(la) 




(lb) 
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study the time-domain behavior of the energy transport in dendritic systems. It can be 
expressed through the total population outside the trap, ip = Yl'iPi^ which is hereafter 
referred to as the survival probability with respect to both decay and trapping. Indeed, 
from Eqs. ([Q) it follows that YliPi = Po + = so that for r{t) one finds 

r{t) = -7^/> -ip . (2) 

Furthermore, the time dependence of r(t) deriving from the decay constant 7, can be ex- 
tracted explicitly by the transformation pi = e~'^^pi. After this transformation Eq. Q is 
reduced to 

r(t) = -e~^'^ = e-^'R{t) , (3) 

where "ifit) = J2[pi, and the pi now obey Eq. (fTb|) with 7 = 0. Thus, \E' is the analog of 
iplt) and represents the total population outside the trap in the absence of excitation decay. 
It is the survival probability with respect to trapping alone. We see from Eq. (jS)) that in 
the time domain, the trapping and the decay of the excitations are independent of each 
other. Therefore, the time behavior of the trapping process can be studied separately, which 
simplifies the analysis. 

The quantity R{t) = — \I' is normalized to unity for finite systems (in sense that 
R(t)dt = 1) and represents the probability distribution of the (pure) trappin g tim e. 
From this, the mean trapping time, often referred to as the mean first passage time is 
calculated in a standard way 

/•oo roo 

{t) = / tR{t)dt = / ^(t)dt . (4) 
Jo Jo 

The inverse quantity {t)~^ represents the effective trapping rate in the absence of excitation 
decay. Note that Eq. (jH) can be rewritten via the Laplace transform R{s) = e~^^R{t) dt 
oiRit): 

I V dR , . 

ds s=o 

which is useful for further considerations (see below). 

We now define the efficiency of trapping (denoted as e) as the total population that is 
transferred to the trap, i.e., the fraction of the initially created excitation that reaches the 
trap during its lifetime: 

/•oo /"OO 

e = lim po{t) = / r(t) dt = / e-^^R{t) dt = R{s)\s=^ . (6) 
Jo Jo 



An other important quantity that contains information about the trapping efficiency is 

r= / ^(t)dt = , (7) 

Jo 7 

which is the mean survival time (with respect to both trapping and decay). Using this 
definition, we can, by convention, define the effective trapping rate W in the presence of 
excitation decay as 

W^ = --7 = T^7, (8) 
r 1 — e 

It should be stressed that e, r, and W , being defined through time-integrations, are 
infiuenced by the excitation decay. In particular, if the latter occurs on a time scale much 
slower than the mean first passage time (t), the trapping efficiency e is close to unity, the 
survival time r is reduced to the mean first passage time (t), given by Eq. (jH), and W ~ {t)^^ ■ 
If however the excitation decay occurs on a time scale that is comparable to or faster than 
the mean first passage time, the quantities e, r, and W will be determined by the interplay 
of the random walk to the trap and the excitation decay. This interplay between trapping 
and decay will be one of the main issues in the remainder of this paper. 

To conclude this section we notice that all quantities relevant to the excitation energy 
transport, such as (t), e, r, and W , are directly related to the Laplace transform R{s) of the 
trapping time distribution R{t) in the absence of decay. In the next section, we will provide 
the exact solution for i?(s). 

III. LAPLACE DOMAIN ANALYSIS: EXACT RESULTS 

From now on, we will consider a specific model for the hopping process, in which only two 
different hopping rates occur. Specifically, we will assume that the hopping rates towards and 
away from the dendrimer's core are, respectively, ki and /c2; no matter at which branching 
point of the dendrimer the excitation resides. This assumption corresponds to the situation 
with a linear energy bias, where the excitation energy difference between units of generation 
M and M — 1 is a constant, AE^ which is identical for every M. As was pointed out in 
Ref. j^, in this case the random walk across the dendrimer may be mapped onto a random 
walk on an asymmetric linear chain, where the rate of hopping is ki in the direction of one 
end of the chain (where the trap resides) and (2; — l)/c2 in the other direction. In other 
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words, instead of a dendrimer of generation a linear chain of length + 1 is considered 
with a trap at site 0. This mapping is illustrated in Fig. |21 

After this mapping, the set of rate equations is different for dendrimers of one, two, and 
N > 2 generations. For a one-generation dendrimer, only one equation occurs: 

A = -hPi . (9) 

For a dendrimer with two or more generations, we have the following equations: 

Pi = - [ki + {z- 1)^2] Pi + hP2 , (10a) 



Pm = {z- l)hPM-i - [h + {z- 1)^2] Pm + hPM+i , 1< M < iV , (10b) 

Pn = {z- l)hPN-l - hPN , (lOc) 

where in case N = 2, Eq. ()10b|) is absent. In all these equations, Pm denotes the total 
population in the Mth dendrimer generation, while the factor z ~1 accounts for the number 
of nearest-neighbor units towards the periphery for each branching point. From this form it is 
clearly seen that the branching {z > 2) leads to a "geometrical" bias towards the dendrimer's 
periphery, even in the absence of an energetic bias {AE = 0). Whether a net bias exists and 
in what direction, depends on the quantity k, = {z — l)k2/ki = {z — 1) exp{—AE /ksT). At 
K, = 1 {ki = {z — 1)^2), the geometrical and energetic biases exactly compensate each other, 
while for /t < 1 (k > 1) a net bias towards (away) from the trap occurs. 

As initial condition we will consider the situation where one excitation has been created 
at the periphery of the dendrimer, i.e., Pa/(0) = Smn- According to Eq. ((Taj), the trapping 
rate in the absence of excitation decay is now given by R{t) = kiPi(t). This will be the 
quantity of our prime interest, for which we will seek a solution in the remainder of this 
section. Solving the one- and two-generation dendrimer problems is straightforward and will 
be done later on. Our main goal is to find the solution of the general problem of an A^- 
generation dendrimer. This may be done in the Laplace domain. If for brevity we introduce 
the dimensionless time t' = kit, Eqs. ()10|1 written in the Laplace domain take the form: 

= -(l + K + s)Pi + P2 , (11a) 
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= kPm^i -{1 + k + s)Pm + Pm+1, km <N , 



(lib) 



-1 = - {I + s) Pn , (11c) 

where the Laplace parameter s is now in units of ki and R = Pi (we use a tilde to denote 
the Laplace transformed functions). Below, we find a recursive relation for R, connecting 
this quantity for dendrimers of different numbers of generations (i.e., lengths of the effective 
linear chain). Therefore, we will from now on denote R for a dendrimer of generations as 
Rn- 

After N — 2 steps of eliminating P/v, Pn-i, -P3 from Eqs. ()llb|l and pic|) . we arrive at 
two coupled equations 

= -{l + K + s)Pi+P2 , (12a) 

-Br, = A^Pi - P2 , (12b) 
a solution of which with respect to Pi = i^^f is 

Rn = , , . (13) 

1 + K + s - An 

Here, A^- and are functions of the Laplace parameter s, which will be specified later on. 
For a dendrimer of + 1 generations, similarly, N — 2 steps of excluding Pn+i, -P/v, -P4 
from Eqs. (jllbjl and ()llc|l yield a system of three coupled equations 

= -{l + K + s)Pi+P2 , (14a) 
= kPi-{1 + k + s)P2 + P3 , (14b) 

-Bn = AnP2 - P3 , (14c) 

where Aj^ and Bjy are the same as in Eqs. ()12a|) . ()12b|) . and (fT^. while now Pi = Rn+i- 
Solving Eqs. ()14a|l - ()14c|) with respect to Pi leads to an expression that is algebraically 
identical to Eq. (fT^ . except that A^ and B^ are replaced by A^+i and B^+i, respectively. 
The latter are given by 

1 + K + s — An 



and represent the recursive relations for these two functions. From comparison of Eq. (jl5|) 
with Eq. (fT^ . one finds that -Bat+i = Rn and A^+i = kR^ /Rn-i- Substituting these 
relations back into Eq. ()13|) . we arrive at a recursive relation for R^: 

J- = (l + K + s)^ . (16) 

Note that Eq. (fTHj) can be used for > 2 only. In order to close the recursive iteration, 
Ri and R2 must be calculated separately. In the Laplace domain the solution of Eq. Q for 

= 1 and Eqs. pOa|) and pOc|) for = 2 is straightforward and gives us the necessary 
quantities: 

Using Eqs. ()17a|) and ()17b|l in the recursive relation Eq. (|T6|l and solving this relation, we 
finally obtain Rn{s) in closed form: 

2^+^ f^^ 1-K + s \ . rr^. 



Rn V V(l + /t + s)2-4 

+ I 1 - — — ^ ~ =] (l + K + s - + K + sy - AkV . (18) 



a/(1 + k + s)2 - Ak^ 

This is the principal result of this section. 

It should be noticed that, despite the fact that Eq. (fTS|) contains square roots, i?^^ 
in reality is a polynomial of A^th order in s. This may be checked by explicit expansion of 
Eq. (fTHj) in powers of s. Alternatively, the recursive relation Eq. (fTBj) itself already represents 
a proof of this statement. 



IV. TRAPPING TIME DISTRIBUTION 



It turns out to be impossible to perform the transformation back to the time domain for 
the general result Eq. (fTHj) . implying that we do not obtain an explicit expression for the 
trapping time distribution. However, the moments of this distribution and its asymptotic 



form at large do allow for an analytical treatment. In this section, we will address these 
two topics. We will also show that the asymptotic form may already be reached for rather 
small dendrimers, making these analytical expressions of practical use. 



A. Moments 

To study the moments of the trapping time distribution, we make use of the following 
relation: 

eR^it)dt = i-ir^ (19) 

Furthermore, we recall that 1/Rn is a polynomial of A^th order in s (see the previous section), 
and thus can be written as 

1 ^ 1 



m=0 ^' 



Substituting Eq. pUjl into Eq. (fTI)|) and comparing coefficients related to the same power of 
s, we obtain a recursive relation for a^-* 

> + {! + (21) 

where it is implied that a}f' = ioi m > N and m < 0. From Eqs. ()17a|l and ()17bj) it 
follows that for = 1 and = 2, 



2 + K . (22) 



(0) (1) (0) (2) , (1) 

al = a\ = = a2 = I , 

This allows us to start the recursive procedure for a!}f' . 

The first important conclusion concerning the expansion coefficients, which follows from 
Eq. (j2H), is that = 1 for any A^. This simply means that the zero'th moment of the 
trapping time distribution -RAr(O) = 1, i.e., the distribution is normalized. As a consequence 
of this fact, one can elucidate the physical meaning of the second coefficient of the expansion 
Eq. fOn|) . It appears to be equal to the first moment of the distribution i?7v, which is nothing 
but the mean trapping time {t)N'- 

1 ARn 



(1) _ A A 



dsR 



N 



s=o \ Rj^ ds 



s=0 



{t)N . (23) 



Substituting this result in Eq. (j^U, yields a recursive formula for {t)^ 

{t)N = 1 + (1 + k) {t)N-l - n{t)N-2 , (24) 
10 



from which earher results for the mean trapping time in dendrimers are reproduced 



(t)iv = ^iV(iV + l) , K = l, (25a) 

N K 

(t)^ = ^ + ^^3^K-l] , n^l. (25b) 

(2) 

The next coefficient, a]^ , contains information about the second moment of R^. It is 



straightforward to show that 



4'' = {t)l-l{t')N. (26) 



On the other hand, using Eq. ()2ip and the method of induction, one may show that 



a!J^ = ^(iV-l)iV(iV + l)(iV + 2), /t = l, (27a) 



(2) _ K^+^[KiN-l)-{N + 2)] 



N{N + 1) n\N + l) + 2K-N 
+ 203^+ (l^^ ' (27b) 



Making now use of Eqs. ^ and ^ in Eq. we find {t'^)N- 

{t^)M = ^N{N + 1)[5N{N + 1) + 2] , K=l, (28a) 

_ NjN + l) iQN + 2)K^+^ 
)^ - (i-«:)2 + (1 - «:)3 

+ -^[«:^^+i + «:^-2], (28b) 
(1 - k)4 

It is worthwhile to consider the asymptotic behavior of the first and second moments of 
Rn at large A^, the case of our primary interest. They are given by 



N , N(N+1) , , 



{t)^ = \N\ (t% = AiV\ «: = !, (29b) 
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These formulas already allow us to make a prediction concerning the shape of the trapping 
time distribution. In particular, comparing Eqs. ()29b|l and ()29c|) with Eq. ()29a|) . we conclude 
that a drastic difference must exist between the two cases k, < 1 (bias towards core) and k > 1 
(no bias or bias away from core), with respect to the shape of the trapping time distribution 
Rn. Indeed, calculating the standard deviation = \/ (t^) n — {t)'^, we obtain 

(Tat = , K < 1 , (30a) 



Ar2 

(Tat = —= , K = 1 , (30b) 



\/6 

= ^ ' fi; > 1 • (30c) 

(1 — kY 

As is seen, for large A^, oj^ <^ {t)^ for k < 1, i.e., the distribution is narrow in the sense 
that its standard deviation is much smaller than its mean. By contrast, ctat ~ {t)^ if k > 1, 
which means that R^ is a broad distribution in this sense. 



B. Asymptotic behavior, N ^ 1 

From Eqs. (j29|) it follows that the characteristic time of trapping {t)]sr in dendrimers of 
higher number of generation is long on the scale of the "effective" hopping times, which are 
(1 — k)"^, {k — 1)^^, and 2 [or {ki — k2)~^, {k2 — ki)~^, and (2ki)~'^ in dimensional units] for 
dendrimers with a total bias towards the trap (k < 1), towards the periphery (k > 1), and 
no bias (k = 1), respectively. Within the Laplace domain, this means that the dominant 
region of the parameter s, being of the order of (t)^^, is, respectively, small compared to 
1 — K, K — l, and 2 for these three different situations. This allows us to significantly simplify 
the expression Eq. (fTH|) for R^. 



1. Total bias towards the trap, k < 1 

We start analyzing the case of a total bias towards the trap (k < 1). Using a Taylor 
expansion of Eq. (fTHj) with respect to s/(l — k), we obtain 

1 f S f S Y 

+1 + ^ • (31) 



R^ [1-kY \ l-Kj V 1 - «^ 
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Since k < 1 and ^ 1, the first term on the left-hand side can be neglected as compared 
to the second one, thus providing us with a very simple expression for 

which can be easily transformed back to the time domain. The result reads 

(1 - I^Y ,N-1 ^-(l-K)i foo\ 

" (iV-1)! • ^ ^ 

We stress that this result is exact in the limit k — (independent of A^). As is seen, R]\f is 
strongly nonexponential and for A^ ^ 1 is characterized by a sharp profile, consistent with 
our findings in the previous subsection. In fact, in the limit of — k) ^ 1 and A^ ^ oo we 
can approximately write = exp[—s{t) jy], where {t)^- = N/{1 — k) is the mean trapping 
time for k < 1. The time-domain behavior, which corresponds to this Laplace transform, 
is Rn = Sit — {t)]y). The nonexponentiality found is in fact a characteristic property of 
the trapping in dendrimers with a total bias towards the trap (k < 1), independent of the 
number of generations A^. 

In order to illustrate how the approximate expression Eq. ()33|) fits the exact result ob- 
tained by numerically integrating Eqs. (fTUj). we plotted in Figs. El and El the distribution R^ 
calculated for k = 1/5 and k = 1/2 at different number of generations A^. Here the solid 
lines are the exact solutions and the dashed lines correspond to the approximation Eq. ()33j] . 
From these plots we conclude that for n = 1/5, Eq. ()33|) works well, even for dendrimers 
with only four generations. On the other hand, for n = 1/2 the deviation of Eq. (jHHI) 
from the exact solution gets larger. These figures also clearly demonstrate the tendency of 
the trapping time distribution for large dendrimers to tend towards a delta-function at the 
mean- first passage time. 



2. Zero total bias, k = 1 

At zero total bias, i.e., when the energetic bias and geometrical one compensate each other 
(k = 1), the problem we are dealing with is equivalent with the classical one-dimensional 
diffusion problem on a finite segment with absorbing and reflecting boundary conditions at 
X = and x = N, respectively, and an initial condition corresponding to the creation of a 
diffusing object at x = A^. The Laplace transform i^^r, derived from Eq. ()18|) in the limit of 
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s <^ 2, reads 

Rn = M— M ~ , ^ , (34) 

+ cosh(v/2i(t)^) ' 

where {t)^- = is the mean trapping time in the diffusive regime of the random walk, 

and we used the fact that (1 ± y/s)^ ~ e^^^ for s <^ 1 and ^ 1. In the time domain 
we then obtain 



2{t) 



N 



(35) 



n=0 

In Fig. El (upper panel), we depicted RN(t) calculated by numertically integrating 
Eqs. for dendrimers of different numbers of generations with k, = 1. First, we note 
that the curves obtained for = 4 and = 8 are almost identical. Second, the dotted 
curve in the plot, corresponding to = 8, coincides in fact with the limiting curve; it is 
not changed by further increasing A^. The decaying part of this curve is nicely fitted by the 
first term of the series Eq. (jHSj), i-e., by exp(— 7r^t/8(t)Ar). Thus, we conclude that already 
for small A^, the approximate result Eq. (jHSl) gives a good fit to the exact result, and the 
tail of RN{t) is described by a single exponential. 



3. Total bias towards the periphery, k > 1 

We proceed similarly to the above in the case of a total bias towards the dendrimer 
periphery (k > 1). Assuming now in Eq. ()18p that s/(fi; — 1) ^ 1, one finds 



We recognize here (t) n = k^^^/ (k — 1)^ as the mean trapping time for k > 1 [see Eq. fl29c|) ]. 
As the relevant region of the Laplace parameter is determined by s(t)Ar ~ 1, while {t)^ ^ 
N/ (/c — 1) at K > 1, the terms s/{k — 1) in the parentheses of Eq. ()36|) are negligible. Upon 
this simplification, Rjsi takes the form 

= Wi^i ' ^^^^ 

which, converted to the time domain, corresponds to the exponential behavior 

1 r M 

~ 77^ exp --— . (38) 

[t/N I \t)N_ 
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Figure El (lower panel) illustrates this finding. All curves presented in this panel were ob- 
tained by numerically integrating Eqs. (fTUI) . As is seen, all curves are close to each other, 
including the one for N = 2. The dotted curve is fitted very well by Eq. (|38|) . From this we 
conclude that the approximation Eq. (j38j) works perfectly for any number of generations. 

The difference in the behavior of (exponential or nonexponential) for different signs 
of the total bias in principle may be used to experimentally probe for the direction of the 
energetic bias in a dendrimer. Indeed, the kinetics of the fluorescence intensity I{t) is 



proportional to ^ = e" 



-7t 



j'Rit')dt' 



i.e., is determined by two decay channels: the 
7 processes and the trapping at the core. Having found the former from, for instance, the 
early-time decay of the intensity, we can then extract information about the direction of the 
energetic bias by measuring I{t). 



V. EFFICIENCY OF TRAPPING 

In this section, we turn to analyzing the trapping efficiency En = RN{s)\s=-y (cf. Eq. Ej), 
the exact expression for which follows directly from Eq. (|T^ : 

= 1+ , ' l + ft: + 7+ V(l + /t + 7)^-4/t 

\ ^(l + «: + 7)2_4«:y V ' J 

+ (l- ../""'M , ) (l + ^ + 7- V(l + ^ + 7)^-4^)'^. (39) 
\ v/(l + K + 7)2-4Ky V / 

From this result, one easily generates plots for any set of variables k, 7, and A^. In partic- 
ular, we present in Fig. IHlthe results for the trapping efficiency according to Eq. (jH^ . for 
dendrimers of A^ = 1 to A^ = 10 generations and different directions of the total bias, setting 
7 = 0.01 (in units of ki). The general trends displayed in this figure are easily understood. 
For small dendrimers, the efficiency is close to unity, unless a strong bias towards the pe- 
riphery is combined with a fast decay. For larger dendrimers, the efficiency decreases due 
to an increased chance of decay before reaching the core. This effect grows when increasing 
the number of generations. However, a more detailed physical interpretation of Eq. (|39|) for 
larger A^ values requires a deeper analysis. In the following, we will focus on this large-A^ 
region. It is then natural to limit ourselves to a decay rate 7 that is small compared to the 
"effective hopping rates" 1 — k, k — I, and 2, as we also assumed with respect to the Laplace 
parameter s in our analysis of the trapping time distribution [see Sec. lIVBj . Otherwise, 
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the trapping efficiency will be low even for a dendrimer of a small number of generations. 
Hereafter, we impose this condition, which allows us to directly use Eqs. fjH^ . and 
(|H7|) . replacing s by 7. 



1. Total bias towards the trap, k < 1 

We ffist consider the case of a total bias towards the trap (k < 1). The corresponding 
expression for is 

= (1 + = exp (-3^) = exp [-7(t)iv] • (40) 

As is seen from this equation, the only parameter that determines the trapping is '^{t)^ = 
'yN/{l — k). The interplay of trapping in the absence of excitation decay and the excitation 
decay itself determines the trapping efficiency: the latter is high (close to unity) for 7(t)Ar <^ 
1, decreasing linearly with 7(t)Ar, and exponentially small in the opposite limit, 7(t)Af ^ 1. 
We note that Eq. implies that in the large- limit with k < 1, (exp[— 7t]) = exp[— 7(t)], 
which is due to the fact that under these conditions the trapping time distribution tends to 
a delta function at the mean trapping time, as we have found in Sec. IIVBI 

Using the definitions Eqs. ((Tj) and (jH)), we can also calculate the mean survival time tj^ 
and the effective trapping rate Wn- They are given by 

1 



tn = - 

7 L 



exp(-7(t)jv) , (41a) 



= ^ . (41b) 

exp [-f{t}N) - 1 

If the decay is slow on the scale of the mean trapping time {'^{t)^ ^ 1), the survival time 
Tn coincides with the mean trapping time {t)]\f, and Wn is just the inverse value {t)]^^. In 
the opposite limit {■y{t)N ^ 1), we get tn = 7~^ (because the excitation cannot reach the 
trap within its lifetime) and Wn = 7exp[— 7(t)7v]- Note that the range of variation of r^r 
is always from {t)^ to 7^^ upon increasing the driving parameter '~f{t)N from zero to values 
large compared to unity. This is the general behavior of the survival time tn, independent 
of the direction of the total bias. 
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2. Zero total bias, k = 1 



In the diffusive regime (k = 1), according to Eq. 



and consequently 



£n = ] , (42) 

cosh y 27(t)Ar 



Wn = , ^, , (43) 

2sinh2v/^(tW2 ' 



where now {t)^ = is the trapping time for k = 1. In the limit of slow decay on the 

time scale of trapping {■y{t)N ^ 1) one obtains 

1 1 

= l + 7(t)^ " l + 7iVV2 ' ^^^^ 
i.e., the trapping efficiency is close to unity, as expected. In Eq. we kept the "small" 
term in the denominator, because this expression works well even if 'j{t)N is slightly larger 
than unity. This is due to the fact that the Taylor expansion of the hyperbolic cosine only 
contains even powers of its argument. 

If '~f{t)N gets larger than unity, another regime of trapping comes into play: 



£iv = 2 exp 



V27(t)^ = exp ( - ^N) . (45) 



As is seen, asymptotically En decreases exponentially, with an exponent proportional to 
(and not A^^). This thus resembles the behavior of trapping in the case of a large bias 
towards the trap (cf. Eq. ()40j). The pre-factor in Eq. ()45|). however, is larger than in the 
case of a bias towards to the trap (7). This is not surprising, because in the diffusive regime, 
the excitation makes steps towards the periphery that slow down the process of reaching the 
trap, allowing for a larger effect of excitation decay before trapping may occur. 

Finally, the effective trapping rate Wn starts from the value {t)jj^ in the limit of slow 
decay {'j{t)N ^ 1) and reveals the same behavior as for 'j{t)]\f ^ 1. 

3. Total bias towards the periphery, k > 1 

In a certain sense, this is the simplest case. For s = 7, Eq. (jTTj) yields for the efficiency 
of trapping 

eN = T— ^77^ (46) 
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for any ratio of {t)^ = k^^^ /{k — 1)^ and 7 ^. Consequently, the relationship Wn = {t)^ 
holds in general. 

VI. CONCLUSIONS 

In this paper we theoretically studied the trapping of excitations in dendritic systems in 
the presence of (radiative and nonradiative) excitation decay when moving to the trap at 
the dendrimer's center. We derived an exact expression for the Laplace transform of the 
trapping time distribution for a dendrimer of any number of generations. This expression 
was then used to analyze the general properties of pure trapping (in the absence of decay), 
focusing on dendrimers of a large number of generations. We found that the general nature 
of this distribution is governed by the total (geometrical and energetic) bias towards the 
trap. In the presence of a bias towards the trap, the trapping time distribution is narrow, 
in the sense that its standard deviation is small compared to its mean. The shape of the 
distribution is strongly nonexponential. Oppositely, in the presence of a bias towards the 
dendrimer's periphery, the trapping time distribution is broad (its standard deviation is of 
the order of its mean), and its shape is essentially exponential with an exponent equal to 
the mean trapping time. The strong difference between both regimes is nicely illustrated by 
comparing Figs. El and El As the fluorescence kinetics is proportional to the trapping time 
distribution, the nature of the fluorescence decay (exponential or nonexponential) may in 
practice be used to distinguish the direction of the energetic bias in dendrimers. We also note 
that, although in the analysis of the trapping time distribution we used the limit of a large 
number of generations, it appears from comparison to numerically exact results that many 
of the analytical expressions hold even for small dendrimers, with just a few generations. 

The trapping efficiency e was found to depend on the ratio of the decay time (7"^) and 
the trapping time in the absence of decay ((t)). The product x = 7(t) turns out to be the 
only essential parameter that governs the trapping in the presence of excitation decay. For 
dendrimers with a total bias towards the trap (k < 1), the efficiency of trapping depends 
exponentially on this parameter within its entire range of values: e = . In the diffusive 
regime of trapping, when the geometrical and energetic bias compensate each other (k = 1), 
e = 1 — X for X < 1, while for x > 1 this behavior changes to a stretched-exponential one, 
e ~ e~^. For dendrimers with a total bias towards the periphery [k > 1), the trapping 
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efficiency e — {1 + x)~^ , independent of x. 

We finally notice that in practice the various regimes with regards to the total bias pa- 
rameter K distinguished by us, may be probed experimentally in one and the same dendritic 
system by varying the temperature. In particular, k tends from zero at low temperature to 
z — 1 aX. high temperature. 
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FIG. 1: Schematic pictures of dendrimers of one (Di), two (D2), and three (D3) generations with 
the coordination number z = 3 (equal to the number of branches at each branching point). The 
large circle at the dendrimer center represents a trap, while the small circles are building units of 
the dendrimer branches, representing sites on which optical excitations can reside. Connected by 
the dotted lines are the dendrimer units belonging to the same generation. 
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FIG. 2: Example of mapping of a dendrimer D3 with z = 3 and a linear energy bias onto 
an equivalent linear chain. The large circle represents the trap, while all units of generation 
M(= 1,2,3) are mapped onto one site M of the linear chain, which is drawn by a small circle. 
The quantities ki and k2 are the hopping rates towards and away from the trap, respectively, in 
the real dendrimer. The factor z—1 counts the number of branches towards the periphery at each 
branching point and multiplies ki to obtain the effective outward hopping rate in the equivalent 
linear chain. 
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FIG. 3: Plots of the trapping time distribution i?Af(t) for dendrimers of different number of 
generations N with a total bias towards the trap (k = 1/5). In all plots, the solid curves represent 
the exact solution for i?Ar(t) obtained from numerically solving Eqs. H10|l . while the dashed curves 
correspond to the approximate expression Eq. H33|). 
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FIG. 4: Same as in Fig. |31 using k = 1/2. 
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FIG. 5: Plots of the trapping time distribution i?iv(t) for dendrimers of different number of 
generations N in the diffusive regime of hopping (k = 1) and for a total bias towards the periphery 
(k = 2). All curves were obtained by numerically solving Eqs. UlUj) . In both panels, the solid curves 
correspond to = 2, the dashed curves to = 4, and the dotted curves to A^ = 8. The decaying 
part of the latter curve in the upper panel coincides with exp[— 7r^t/8(t)Af]. 
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FIG. 6: Plots of the trapping efficiency ejv as a function of number of generations A'" calculated 
for 7 = 0.01 and different values for k. The squares correspond to k = 1/2 (total bias towards the 
trap), circles correspond to k = 1 (no total bias, diffusive regime of hopping), and the triangles 
correspond to k = 2 (total bias towards the periphery). 
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